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Abstract. We develop direct and inverse spectral analysis for fi- 
nite and semi-infinite non-self-adjoint Jacobi matrices with a rank 
one imaginary part. It is shown that given a set of n not neces- 
sarily distinct non-real numbers in the open upper (lower) half- 
plane uniquely determines a n x n Jacobi matrix with a rank one 
imaginary part having those numbers as its eigenvalues counting 
multiplicity. An algorithm for reconstruction for such finite Jacobi 
matrices is presented. A new model complementing the well known 
Livsic triangular model for bounded linear operators with rank one 
imaginary part is obtained. It turns out that the model operator 
is a non-self- adjoint Jacobi matrix and it follows from the fact that 
any bounded, prime, non-self-adjoint linear operator with rank one 
imaginary part acting on some finite-dimensional (resp., separable 
infinite-dimensional Hilbert space) is unitary equivalent to a finite 
(resp., semi-infinite) non-self-adjoint Jacobi matrix. This result 
strengthens the classical Stone theorem established for self-adjoint 
operators with simple spectrum. We establish the non-self-adjoint 
analogs of the Hochstadt and Gesztesy-Simon uniqueness theo- 
rems for finite Jacobi matrices with non-real eigenvalues as well as 
an extension and refinement of these theorems for finite non-self- 
adjoint tri-diagonal matrices to the case of mixed eigenvalues, real 
and non-real. A unique Jacobi matrix, unitarily equivalent to the 
operator of integration (Ff)(x) = 2 i f f(t)dt in the Hilbert space 
C%[0, 1] is found as well as spectral properties of its perturbations 
and connections with well known Bernoulli numbers. We also give 
the analytic characterization of the Weyl functions of dissipative 
Jacobi matrices with a rank one imaginary part. 
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1. Introduction 



Self-adjoint (or real) finite and semi- infinite Jacobi matrices of the 
form 
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where > 0, and are real numbers for all k = 1,2,... play an 
important role in various problems of mathematical analysis and theo- 
retical and mathematical physics. They appear as the discrete analog 
of Sturm-Liouville operators, in inverse spectral theory, in the study of 
the classical moment problem, and in the investigation of completely 
integrable nonlinear lattices [I], jH], 0, [J2|, [TB] , 
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More general tri-diagonal matrices with complex entries (or complex 
Jacobi matrices) also attracted much attention as a useful tool in the 
study of orthogonal polynomials, in the theory of continued fractions, 
and in numerical analysis 0, [Zj, [HI]- 

Let the linear space C n of columns be equipped by the usual inner 
product 

n 

{ X :V) = ^2 x kVk 
k=l 

and let ^i^) be the Hilbert space of square summable complex- valued 
sequences 

x = {xi,x 2 , ...,x k ,...} 

which we consider as semi-infinite vector-columns with the inner prod- 
uct given by 

oo 
k=l 

Let 

5 k :=(0 1 _^0 1 l,0,...) T } k = l,2,.... 

k 

Then the vectors {S k } form an orthonormal basis in C™ (resp.^N)). 
An n x n complex Jacobi matrix J determines a linear operator in 
the Hilbert space C n by means of the matrix product J ■ x. For the 
semi-infinite case we will suppose in addition that 

(1.3) sup{|a fc | + \h\} < oo, k G N. 

k 

This condition is necessary and sufficient [2] for boundedness of the 
Jacobi operator in ^(N) defined as J ■ x, where J is a semi-infinite 
complex Jacobi matrix. Moreover, under the condition a k ^ for all 
k, the vector 5± is cyclic for the Jacobi operator J. The classical Stone 
theorem [T], j2], [22] states that every self-adjoint operator with simple 
spectrum in a separable Hilbert space TL is unitarily equivalent to the 
operator determined by a self-adjoint Jacobi matrix. 

In this paper we consider finite and semi-infinite tri-diagonal matrices 
of the form (IT~T|) and JOJ with 

(1.4) Imfox > 0, b k = b k (k = 2,3,...), a k > (k G N). 

Such matrices determine bounded Jacobi operators J in C n or in ^(N) 
that posses the properties 

Im (Jx, x) > for all x G C n or x G l 2 (N), 

ran (J - J*) = {\8 l , X G C}. 
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In the following we will call such tri-diagonal matrices the dissipative 
Jacobi matrices with rank one imaginary part. Because the vector 8\ 
is cyclic for J, the dissipative operator J is prime [TUj, [H], [21!, and 
therefore it has no real eigenvalues. 

M. S. Livsic |2Zj (see also jTTj) constructed a triangular model for 
a bounded prime dissipative operator A with a rank one imaginary 
part. The method of construction is based on the factorization of the 
characteristic function W(z) of the operator A. The model operator is 
in general a coupling of two triangular operators. The first one is given 
by a finite or semi-infinite triangular matrix of the form 



A 



a 2 + 



a 



n + 2 PnJ 



or 



A 



(a x + §/?? 1(3x132 
a 2 + \i3\ 



\ 



iPiPk 
if32f3k 



where {ak} k =i are real numbers and {/3fc}'£ =1 , are positive numbers and 
Zk — a k + i(P%/2), k = 1, 2, ... are the non-real eigenvalues of A. The 
second one is the operator in the Hilbert space C 2 [0, 1} given by 



(Bf)(x) = a(x)f(x)+i / f(t)dt 



where a is a bounded nondecreasing right continuous function on [0, 1]. 
The operator B has pure real spectrum which coincides with the closure 
of the set {a(x), x G [0,/]}. 

One of the main results of the paper states that a bounded prime 
dissipative operator with a rank one imaginary part acting on a sepa- 
rable Hilbert space is unitarily equivalent to the operator determined by 
a dissipative Jacobi matrix with rank one imaginary part. This result 
is a non-self-adjoint analog of the classical Stone theorem, established 
for self-adjoint operators with simple spectrum in pQ 7 [2], |H2]- Thus, 
dissipative Jacobi matrices with rank one imaginary parts provide new 
models for the prime bounded linear operators with rank one imaginary 
parts. The entries of the corresponding Jacobi matrix can be found 
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using the continued fraction (J-fraction) expansion [34"j 



M(z) 



z-b 1 + z-b 2 + z~b 3 +... + z — b n + ... ' 



where W(z) is the characteristic function, M(z) = ^(W(z) — 1), and 
P = lim (iz{l - W(z))) . 



By the Livsic theorem |2Z], a prime dissipative Volterra operator A 
acting on the Hilbert space H, with rank one imaginary part such that 



for all h & Tl, ||e|| = 1, is unitarily equivalent to the integration oper- 
ator of the form 



acting on the Hilbert space /^[O, I]- In this paper a unique Jacobi matrix 
unitarily equivalent to the integration operator T is found, and some 
connections with the well known Bernoulli numbers are established. 

We also show that for n not necessarily distinct given complex num- 
bers in the open upper half-plane, there exists a unique n x n Jacobi 
matrix having those numbers as its eigenvalues counting multiplicity. 
Two algorithms for the reconstruction of a finite Jacobi matrix of the 
form from its eigenvalues are presented. We obtain the following 
non- self- adjoint analogs of the Hochstadt and Gesztesy-Simon unique- 
ness theorems for self-adjoint finite Jacobi matrices: Let J be an n x n 
Jacobi matrix satisfying conditions ([1.4)1 . Suppose that the eigenval- 
ues zi, . . . , Zk (taken from the open upper half-plane) with their mul- 
tiplicities h, . . . , Ik are known together with bi, a%, b 2 , . . . , a n _ r _i, b n - r , 
where r = li + Z 2 + • • • + h- Then a n _ r , 6„_ r+ i, . . . , a„_i, b n are uniquely 
determined. This as well as some extensions and refinements of the 
Hochstadt and Gesztesy-Simon uniqueness theorems concerning the case 
of dissipative but non prime tri-diagonal matrices with rank one imag- 
inary part are established. 

Our main tools are the Livsic characteristic function, its linear- 
fractional transformation jTUj, jTTj . |27| and the Weyl function of 
a Jacobi matrix. Moreover, in order to establish the above mentioned 
non-self-adjoint analogs of the Hochstadt ^J] and Gesztesy-Simon jT3] 
uniqueness theorems in the theory of inverse spectral problems we de- 
velop further the approach Gesztesy-Simon used in [To] . 
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(A-A*)h = 2il(h,e)e 
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Our paper is organized as follows: In Section|2]we present some prop- 
erties of prime, bounded, dissipative operators with a rank one imag- 
inary part and the Livsic theory of their characteristic functions and 
triangular models. We define the Weyl function of such operators and 
give its analytic characterizations. We emphasize that the important 
role of the Weyl function (Weyl-Titchmarsh function) in the spectral 
theory of self-adjoint differential and difference operators and its appli- 
cations to nonlinear equations is well-known (see |T], [5], |H1, M [T5] . 
[To] . [33]). In Section El some basic properties of complex Jacobi matri- 
ces and the corresponding Jacobi operators as well as a survey of the 
inverse spectral problems for finite and semi-infinite self-adjoint Jacobi 
matrices are given. We begin our study of dissipative Jacobi matrices 
with a rank one imaginary part in Section 0] In particular, the connec- 
tion between m+- and m_-functions of a finite dissipative Jacobi matrix 
is established. In Section El we present solutions of the inverse spectral 
problems for finite dissipative Jacobi matrices. The non-self-adjoint 
analogs of the Hochstadt and Gesztesy-Simon theorems are presented 
in Subsections 15.31 and 15.41 The non-self-adjoint analog of Stone's the- 
orem is considered in Section El We also obtain in Section^ that if the 
Livsic characteristic function W(z) of a bounded dissipative operator 
with rank one imaginary part satisfies the condition W(— z) = W~ 1 (z) 
in some neighborhood of infinity, then the corresponding Jacobi matrix 
possesses the property Re&i = &2 = • • • = 0. In Section^ we find the 
Jacobi matrix corresponding to the integration operator and establish 
that the perturbation of its upper left entry b\ leads to a Jacobi matrix 
with a complete system of eigensubspaces. 

We will use the following notation: For a bounded linear operator A 
in 7i, Ar and Aj denote its Hermitian components 

A + A* „ A - A* 
Ar = — - — , Aj = 



2i 1 

and p(A) and cr(A) denote the resolvent set and the spectrum of A. By 
C + (resp.,C_) we denote the open upper (resp., lower) half-plane. For 
a continued fraction 

Cl 

dl + — C2 



do 



d 3 + • • ■ 

we use the notation (see [23], |S1- [21]) 

c i C2 C3 

di + d 2 + d 3 + . . . 
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Finally, we would like to mention that the fundamental research pro- 
vided by Yuri Berezanskii, Harry Hochstadt, Fritz Gesztesy and Barry 
Simon in the theory of self-adjoint Jacobi matrices as well as by Moshe 
Livsic in the theory of non- self- adjoint operators inspired and encour- 
aged us to make a new step and to develop the direct and inverse 
spectral analysis of Gnite and semi-infinite non- self- adjoint Jacobi and 
tri-diagonal matrices with a rank one imaginary part as presented in 
this paper. 



2. Bounded prime dissipative operators and their 
characteristic functions 

2.1. Bounded prime dissipative operators. Let 7i be a complex 
separable Hilbert space with the inner product (•, •) and the norm 1 1 • 1 1. 
The operator A is called dissipative if A% > 0. The resolvent set p(A) 
of dissipative operator A contains the open lower half-plane and the 
estimate 



WiA-ziy^l < -r^—r, Imz < 



holds. Let 



(2.1) H s = span {A n A I U,n = 0,1,...}. 

The subspace Ti s reduces A and A \ (Ti.QT~C s ) is a self-adjoint operator. 
The operator A\7i s is called the prime part of A. 

Definition 2.1. A bounded linear operator A in a separable Hilbert 
space is called prime or completely non-self-adjoint if there is no re- 
ducing invariant subspace on which the operator A is self-adjoint. 

It is well known JU] that an operator A is prime in TC if and only if 

(2.2) spm{A n A I H,n = 0,l,...} = H. 

Proposition 2.2. Let A be a dissipative operator with a rank one imag- 
inary part and let g be a vector in T~C such that 

(2.3) 2A//i = (h,g)g, heH. 

Then A is prime if and only if the vector g is cyclic for the real part 
A R . 

Proof. Suppose that A is prime. Then ()2.2|) holds. Let us prove that g 
is a cyclic vector for Ar. Let 

H' = span {{A R ) n g, n = 0, 1, . . .} ^ H. 
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Then TC and TC" = Tt Tt' are invariant with respect to A R . Since 
Tt" C KerA/, it follows that 

A R \H" = A\H" = A*\TC", 

A* n H" C Tt and A I A* n \H" = for all n = 0, 1, ... . Now from (1231) 
we obtain that W = {0}, i.e. g is a cyclic vector for Ar. 

Conversely, suppose that the vector g is cyclic for Ar, i.e. TC = TC. 
Let the subspace TC S be defined by (|2.1jl . Then 

and Tt QTt s as well as W s reduces Ar. Because g E TC S , we get that 
W C 7i s . It follows that 7i s = TC, i.e. A is a prime operator. □ 

A prime dissipative operator has no real eigenvalues and its non- 
real eigenvalues belong to the open upper half-plane. It is known (see 
|lUj ) that the non-real spectrum of an operator with compact imaginary 
part consists of eigenvalues of finite algebraic multiplicities (dimensions 
of the corresponding root subspaces) and the limit points of non-real 
spectrum belong to the spectrum of the real part of operator. Note that 
any eigensubspace of an operator with a rank one imaginary part which 
corresponds to a non-real eigenvalue has a dimension equal to one. 
Observe also that if dimH = n and A is a dissipative operator with a 
rank one imaginary part Aj, and if (J2.3)) holds, then as it follows from 
Definition 12. 21 and Theorem 12 . 61 the following conditions are equivalent: 

(1) A is a prime operator; 

(2) A has no real eigenvalues; 

(3) the vectors g, Ag, . . . , A n ~ 1 g are linearly independent. 

The following completeness criterion has been established by M. S. Livsic 

m\ 

Theorem 2.3. [27 . Let A be a bounded prime dissipative operator in 
a separable Hilbert space Tt with the finite trace SpAj. Then the closure 
of the linear span of all root subspaces coincides with TC if and only if 

(2.4) ^Imz^Sp^j, 

n 

where {z n } is the set of all non-real eigenvalues of A counting multi- 
plicity. 

2.2. The Livsic characteristic function of non-self-adjoint bounded 
operators. In this subsection we present basic facts of the Livsic char- 
acteristic functions theory jTU], [H], |2ZI, |2Hj- We restrict ourselves 
just to the case of a bounded dissipative operator A with a rank one 
imaginary part Aj. 
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A vector g satisfying ()2.3|) is called a channel vector for A. The 
function 

(2.5) W(z) = 1 - i ((A - zl)- l g, g) , z e 

is called a characteristic function of A. Because an operator A is 
bounded and dissipative the characteristic function W(z) possesses the 
following properties: 

(1) \W(z)\ > 1 for z E p(A), Imz > 0, 

(2) \W(z)\ < 1, Imz < 0, 

(3) |W(z)| = 1 for z e p(A), Imz = 0, 

(4) jy (z) = l+i Xl^Lo ^n+f^ in a neighborhood of infinity. 
In addition the function 



1 _ W(z)W(£) 

(2.6) G(z.,£) = =r is a nonnegative kernel, 



i.e. 



n 



/ J XiXjG(zj, Zi) > 0, (n = 1, . . . ) 

for any set of points {^}" =1 C p(^4) and any vector (Aj, ■ • ■ , A n ) G C n . 

Theorem 2.4. |27j . Lei Ai and &e bounded dissipative operators 
with a rank one imaginary part. If the corresponding characteristic 
functions W\(z) and W^{z) coincide in some neighborhood of infinity, 
then the prime parts of A\ and A 2 are unitarily equivalent. 

Theorem 2.5. [Hi]- Let a function W(z) be holomorphic outside some 
closed bounded subset V of the closed upper half-plane and possesses the 
properties 

(i) \W(z)\ < 1, Imz < 0, 

(ii) \W(z)\ = 1 for z e C\V, Imz = 0, 
(hi) lim W(z) = 1. 

z— *oo 

Then there exists a prime bounded dissipative operator with a rank one 
imaginary part whose characteristic function coincides with W(z). 

Theorem 2.6. 27 . The characteristic function W(z) of a bounded 
dissipative operator A with a rank one imaginary part admits the fol- 
lowing multiplicative representation: 

N ' - \n k ( I 1 

Zk\ I dx 



/ — \ n k I r 



a(x) — z 
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where {zk}^ =1 (N < oo) are the distinct non-real eigenvalues of A, 
{n^} are their algebraic multiplicities, a(x) is a bounded nondecreasing 
and continuous from the right function on the interval [0,1]. 

Observe, that for any operator A and his channel vector g (12. 3 j) 

k=l 

Moreover, from Theorems 12.31 and 12.61 it follows that the closure of the 
linear span of all root subspaces of operator A coincides with TC if and 
only if 

n \\g\\ 2 
k=i 

2.3. Triangular models for bounded dissipative operators with 
rank one imaginary parts. The next theorem provides a triangu- 
lar model for a finite-dimensional prime dissipative operator with one 
dimensional imaginary part. 

Theorem 2.7. [22j. Let dimTC = n and let A be a dissipative 

operator in TC with a rank one imaginary part Aj. Then there exists 
an orthonormal basis in TC in which the matrix of the operator A takes 
the form 



A 



(a x + {pi ifch ■ ■ ■ ifcp n \ 

a 2 + l(3 2 2 ■ ■ ■ il3 2 (3 n 



\ ■ ■ ■ a n + iP 2 J 

where {ak} k= i ore real numbers and {(3k}k=i are nonnegative numbers. 
The operator A is prime if and only if the numbers {(3^} are positive. 
The characteristic function of the operator A takes the form 

w(z) = TT z ~ at + 1131,2 

Now consider the case of infinite-dimensional operator. 

1) Suppose that a bounded prime dissipative operator A with rank 
one imaginary part has a complete system of root subspaces. Then the 
characteristic function of A takes the form 

oo _ 
Z — Zk 



n 



where {z k } are the set of not necessarily distinct complex numbers 
counting multiplicity and such that 
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(1) lmz k > 0, 

( 2 ) Efcli Im ^ < °°> 

(3) \z k \<C. 

From the Livsic approach [27| such an operator A is unitarily equivalent 
to the operator A in the Hilbert space / 2 (N) given by the triangular 
matrix 

/c*i + \ (31 i/3i/3 2 

1^2 








a 2 + m 







V 



where {«fc}^ =1 are real numbers and {Pk}k=i are positive numbers. 

2) Suppose that A has only real spectrum. Then the characteristic 
function of A takes the form 



W(z) = exp | —i J 



dx 



a(x) — z 



Consider the following operator B in the Hilbert space L 2 [0,/]: 

/ 

(Bf)(x)=a(x)f(x)+i J f(t)dt. 

X 

The operator A is unitarily equivalent to the prime part of the operator 

B |2ZI- 

3) Let A be an arbitrary bounded prime dissipative operator with a 
rank one imaginary part acting on infinite-dimensional Hilbert space 7i. 
Then its characteristic function admits multiplicative representation 
(12. 7j) . Consider Hilbert space fj = C n © L 2 [0, 1} in case when the non- 
real spectrum of A is finite and consists of n complex numbers counting 
their algebraic multiplicity and the Hilbert space S) = / 2 (N) © L 2 [0, 1] 
if the non-real spectrum of A is infinite. Let G be an operator in Sj 
given by the block operator matrix 



G 



a r 
o b) ' 



where V : L 2 [0, 1] — > C n or T : L 2 [0, Z] — » / 2 (N) is given by the formula 

/ i 

(r/)(x)=i ; " 
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where 8\ G C n or 8\ G ^(N) is the vector with the first component 
equal to 1 and remaining components equal to zero. Then A is unitarily 
equivalent to the prime part of G fI7\ . 

2.4. Linear-fractional transformation of the characteristic func- 
tion. Let T be a bounded linear operator acting on some Hilbert space 
TC and let h be a nonzero vector in TC. Consider the family of rank one 
perturbations of T given by 

T t = T + t(-,h)h, 

where t is a complex number. One can easily derive the following 
relation 

(2.8) ((T, - h) = • « * PPi) n P(T). 

Moreover, the number G p(T) is the eigenvalue of T t if and only if 
Zq satisfies the equation 

1+t ((T — z I)^ 1 h, h) = 0. 

Let A be a dissipative operator in TC with a rank one imaginary part 
Aj and let g G TC such that (|2.3j) holds. Define the following function 

(2.9) V(z) = \ ((A R - ziy'g, g), ze p{A R ). 

The function V(z) is holomorphic on the domain C \ [a,b], where a 
and b are the lower and upper bounds of the spectrum of self-adjoint 
operator A R , and posseses the properties 

(1) v w:yw>o, zec+ uc_, 

(2) V'(l) = V(z), 

(3) V(z) = —\ ^2 — — ^n+i ^ * n a neighborhood of infinity. 

n=0 Z n 

The properties (1) and (2) mean that V(z) is a Herglotz-Nevanlinna 
function and admits the integral representation 

b dn(t) 

where Q(t) = (E(t)g,g)/2 and E(t) is a resolution of identity for the 
operator A R . Since 

A R = A- i -{-,g)g, 
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from (Q and (J22J) it follows that (HO]) 

(2.10) ^j.zep^n^. 

r l J i+iv(z) 

The non-real spectrum of A consists of all z, Im z > which are solu- 
tions of the equation 

V{z) = i. 

2.5. The Weyl function of a bounded dissipative operator. In 

this subsection we define and study the Weyl function of a prime dissi- 
pative operator with rank one imaginary part. In the classical case of 
the Weyl functions of self-adjoint differential and difference operators 
we refer to PQ, 0, 0, 0, [E], [SSI- 
Definition 2.8. Let A be a prime dissipative operator with a rank one 
imaginary part. Let e G ran (A — A*), \\e\ \ = 1. The functions 

(2.11) m A (z) = ((A - zT)~ x e, e) , z G p(A) 
and 

(2.12) m AR (z) = ((A R - ziy'e, e) , z e p{A R ) 

are said to be the Weyl functions of the operator A and A R , correspond- 
ingly. 

Since e G ran (A — A*) and ||e|| = 1, the operator Aj takes the form 

Ajf = Z(/,e)e, 

where I > 0. Observe that by Proposition 12.21 the vector e is cyclic for 
A and A R . The function m A (z) is holomorphic on the resolvent set of 
A which includes the open lower half-plane and the neighborhood of 
infinity of the form \z\ > \ \A\ \ and in this neighborhood 

(A k e,e) 



z k+i 



k=0 

The function m AR (z) is a Herglotz-Nevanlinna function holomorphic 
on the resolvent set of A R and for \z\ > \\A R \\ 



fc=0 

Relations (Q and (JUJ) yield 

m AR {z) = - V(z), z G p(A R ) 
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From (J2.3j) and (|2.5|l it follows that the characteristic function W(z) of 
A and the Weyl function m A (z) are connected via relations 

(2.13) W{z) = 1 - 2ilm A {z), m A {z) = ^(W(z) - 1), z G p(A). 

The relation ()2.8)1 gives the following connections between the Weyl 
functions of A and Ar 

mA (z) = . ^ . , 

(2.14) 

. . im A (z) , AS / \ \ 

z + lm A [z) 

From (j2.14|) it follows that the non-real eigenvalues of A coincide with 
poles of the function m A (z) and zeroes of the function i — lm AR (z). 
The following theorem directly follows from Theorems 12 .41 and relations 

Theorem 2.9. Lei Ax and fre two prime dissipative operators with 
rank one imaginary parts. If the Weyl functions of A\ and A 2 coincide 
in some neighborhood of infinity, then the operators A\ and A 2 are 
unitarily equivalent. 

Suppose that a Hilbert space 7i is decomposed as Ti = Hi © Ti.2- 
Then every bounded operator A in Ti has the block-operator matrix 
representation with respect to this decomposition 



.4 



An A u 
A21 A 2 2 



According to the Schur-Frobenius formula the relation 

(2.15) P m {A - zI)- l \Hi = (-zl + An - Ai 2 (A 2 2 - zI)~ 1 A 21 )~ 1 

holds for z G p(A) fl p(A 2 2). Here P-^ 1 is the orthogonal projection in 
Ti onto Tii. Let A be a dissipative operator with a rank one imaginary 
part and let Tii = ran (A — A*). Then H.2 = Ker (A — A*), the operator 
A22 is self-adjoint in H2, and A12 = A 21 . It follows that 

A n e = be, A 12 f = a(f, h)e, f G H 2 , A 2i e = ah, 

where e G ran (A — A*) and h G Ker (A — A*) are the unit vectors and 
Imb > 0. Hence (I2.15|) takes the form for z G p(A) fl p{A 2 2) 

(2.16) {(A - zl)~ l e, e) ' 



-z + b- \a\mA22- ziy^^h)' 
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If A is a prime operator, then a^O and in this case for z G p(A)r\p(A 22 ) 

(2,7) ((A 22 - ,1)-% h) = ^ + 6 - {{A _ z ) ) . le>e) ) ■ 
It follows that 

((A 22 - ^7)-^, fc) = Jjj (ft - * - —L^ , 2 g p(A) n p(A 22 ), 

((A 22 - ziy l h, h) = —^[Reb-z ^— ) , z G n p(A 22 ). 

\a\ V m A R {Z) J 

In addition, the number zq, Im^o > is the eigenvalue of A if and only 
if z satisfies the equation 

b-z- \a\ 2 ((A 22 - zI)-% h) = 0. 

The corresponding eigenspace is 

span {e — a(A 22 — zqI)^ 1 ^ . 

The next theorem gives the analytical characterization of the Weyl 
functions of prime dissipative operators with a rank one imaginary 
part. 

Theorem 2.10. Let M(z) be a function holomorphic outside some 
closed bounded subset V of the closed upper half-plane. The following 
statements are equivalent: 

(i) The function M(z) is the Weyl function of some bounded prime 
dissipative operator with a rank one imaginary part. 

(ii) (a) The function M \z) has the asymptotic expansion at infinity 

(2.18) M {z) = ---- + o' ' 



z z 2 \ z 3 



where Imb > 0. 
(b) The function 

W(z) = 1 - 2ilmbM(z) 

is the characteristic function of some dissipative bounded 
operator with a rank one imaginary part. 
(iii) (a) The function M(z) has the expansion ()2.18j) at infinity. 
(b) The function 



(2.19) K( . fl;= M (z )-M(fl + aim t M W M(fl 
is a nonnegative kernel. 
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(iv) (a) The function M(z) has the expansion (|2.18|) at infinity. 
(b) The function 

Q(z) = tM{z) 
^ y> ' i + lmbM{z) 

has analytic continuation onto the outside of some bounded 
interval of the real axis as the Herglotz-Nevanlinna func- 
tion . 

(v) (a) The function M(z) has the expansion (|2.18|) at infinity. 
(b) The function 

has analytic continuation onto the outside of some bounded 
interval of the real axis as the Herglotz-Nevanlinna func- 
tion. 

Proof. Let M{z) be the Weyl function of some bounded prime dissi- 



pative operator A with a rank one imaginary part. From (j2.11|) we 

get 

M(z) = ((A-ziy^e)) = -J2' 1 



k=0 



z k+l 



z z 2 \ z 3 



O — , \z\> \\A\ 



where e G ran (A — A*), \\e\\ — 1, b — (Ae, e), and Imb > 0. So, 
(I2.18j) holds. Since Aie = Imbe, from (|2.5|) it follows that the function 
W(z) = 1 — 2ilmb M(z) coincides with the characteristic function of 
A. Hence 



_ M(z) - M(fl + lmbM(z)M(Q _ 1 - W(z)W(Q 
is a nonnegative kernel. By ()2.14j) the function 

Q{z) = -* M{ h Z L v z G M n Mr) 

1 + imo M(z) 

coincides with the Weyl function of Ar and therefore, has analytic 
continuation onto the outside of some bounded interval of the real axis 
as the Herglotz-Nevanlinna function. By (J2.17|) . the function 
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also has analytic continuation onto the outside of some bounded in- 
terval of the real axis as the Herglotz-Nevanlinna function. Thus, the 
statement (i) implies the statements (ii), (iii), (iv), and (v). 

Now suppose that the statement (ii) holds true, i.e. there exist a 
Hilbert space Ti and a bounded prime dissipative operator A in Tt 
with a rank one imaginary part such that 

1 - i((A - zl)~ l g, g) = l- 2z'Im bM(z), 

where g G ran Aj and Aj = (-,g)g. It follows that \\g\\ 2 = 21mb and 

M(z) — ((A - zl)~ x e, e), where r- fJ 



V2lmb 

By definition 12.81 the function M(z) is the Weyl function of A. So, 
(ii)=> (i). If (iii) holds, then for W{z) = 1 - 2i\mbM{z) from (l2~T3j) 

we obtain that 

1 - W(z)W{£) 

Since /C(z, £) is a nonnegative kernel, we get that |W(2)| < 1, Imz < 0, 
and |W(z)| = 1 for 2 G C \ D, I1112 = 0. In addition, from pT8|) 
we get that lim W(z) = 1. By Theorem 12.51 the function W(z) is 

z— >oc 

the characteristic function of some bounded prime dissipative operator 
with rank one imaginary part. Thus (iii)=^(ii). 
Let us show that (iv)=^(ii). Define the function 

= l-ilmbQ(z) 
{ } l + i\mbQ{z)' 

Then W(z) = 1 — 2ilmbM(z). Since Q(z) is a Herglotz-Nevanlinna 
function, we have \W(z)\ < 1, Imz < 0, |W(z)| = 1 for z e C \ V, 

Imz = 0, and lim W(z) = 1. By Theorem 12.51 the function W(z) is 

2^00 

the characteristic function of some bounded prime dissipative operator 
with a rank one imaginary part. 

Suppose that (v) holds true. Because 

is a Herglotz-Nevanlinna function, the function 



z-l 

is a nonnegative kernel. From ()2.19|) it follows that 



m(z) - m(0 



K{z,t) = M{z)M{i) 



18 



YURY ARLINSKII AND EDUARD TSEKANOVSKII 



Therefore, /C(z, £) is a nonnegative kernel and (v)=>- (ii) 



□ 



3. Jacobi matrices and their Weyl functions 

3.1. Complex Jacobi matrices and corresponding Jacobi op- 
erators. Let {bk}, {at}, cik 7^ 0, k > 1 be complex numbers. A 
tri-diagonal matrix of the form 



(3.1) 



J 



(b x a x 
ai b 2 a 2 
a 2 h a 3 



V 



o 



a>n-i b n / 



is called finite complex Jacobi matrix and a matrix 

(b x ai ■ • -\ 

a x b 2 a 2 • • 
a 2 b 3 a 3 ■ • 



(3.2) 



J 



V 



7 



is called semi-infinite complex Jacobi matrix. The case of real entries 
{bk} and positive entries {a^} corresponds to the classical symmetric 
Jacobi matrix pp, [23- We will call such matrix a self-adjoint 

Jacobi matrix. A finite nxn Jacobi matrix determines a linear operator 
(the Jacobi operator) in the Hilbert space C n . Let Co be the linear 
manifold of vectors in l 2 (N) with finite support. A semi-infinite Jacobi 
matrix determines two linear operators in l 2 (N) given by the formal 
matrix product J -x. The first operator is defined in C - This operator 
is densely defined and is closable. Let [ J] m j n be its closure. The second 
operator [J] maa; has the domain 

dom ([J] max ) = {x E Z 2 (N) : J ■ x E Z 2 (N)}. 

A semi- infinite Jacobi matrix is called proper if [J] m i n = [J] ma x (cf. 
[H]). It is well known (cf. |2], |H]) that dom [J] m i n = l 2 (N) and [J] m i n is 
bounded in l 2 (N) if and only if the entries {a^} and are uniformly 
bounded, i.e. condition (|1.3|) is fulfilled. The tri-diagonal matrix is 
compact if and only if 



lim bi 



lim afc 

fc^oo 



0. 



A Jacobi matrix is called bounded if (jl.3|) holds. 

Because 

{J k di) k+1 = ak^k-i ■■■ax, 
(J%) =0, m> k + 2, 



(3.3) 
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and a k 7^ 0, the vectors Si, JSi, . . . , J k Si, . . . are linearly independent 
and moreover the vector Si is cyclic for the operator [J] m j n in ^(N). 
The system of second order difference equations 



(3.4) 



a k P k+ i(z) + b k P k (z) + a k -iPk-i(z) = zP k (z), k>l,a := 1, 
P o (z)=0, P 1 (z) = l 



determines polynomials P k (z) of degree k — 1 (in the case of n x n 
Jacobi matrix we define a n := 1). Note that 



(3.5) P k+ i{z) 



ai ■ ■ ■ a k 



-z k + 



lower degree in z, k 



n. 



3.2. Finite complex Jacobi matrices. The next proposition prob- 
ably is well-known. We give a proof for completeness. 

Proposition 3.1. Let J be annxn complex Jacobi matrix with a k ^ 
for all k and let polynomial {P k (z)}Q +1 be defined by ()3.4|) . Then z 
is an eigenvalue of J if and only if P n+ i(z ) = 0. Moreover, if zo is a 
root of the polynomial P n+ i(z) of the multiplicity I, then vectors in C n 
defined as follows 



( 1 \ 

^o) 



\Pn(Zo)J 



ei 



(3.6) 



/ 



ei-i 



( o \ 

P^zz) 



\P'n(Zo)J 

o \ 






e 2 



1 
2! 






P^(zo) 



\P^z )J 



pr\zo) 



\Pt 1] (z )J 
satisfy the relations 

(J - z I)e = 0, (J - z I)e k+1 = e k , k = 0, 1, 2. 
In addition, the relation 

(3.7) P k +i{z) =a 1 a 2 ■ ■ ■ a k det (zl - J [ljk] ) , k = 1, 
holds, where J\i jk ] is the k x k upper left corner of J. 



n 
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Proof. Let z be an eigenvalue of J and let (jji y 2 ... y n ) be the 
corresponding eigenvector. Then we have the system of linear equations 

' feii/i + a x y 2 = z y 1 

aiyi + 622/2 + 022/3 = z y 2 
< 

On-22/n-2 + K-lVn-l + Ctn-lVn = ZoVn-1 

k a n _i?/ n _i + b n y n = z y n - 

Since 7^ for all = 1, 2 . . . , n — 1, we can express y 2 > 2/3, • • •> 2/n 
linearly through y x . It follows that 1/1 ^ and we may put y\ = 1. 
Comparing with (|3.4j) we obtain that y 2 = -^2(^0)? • • • > Un = Pn{ z o), 
and P n+ i(z ) = 0. 

Conversely, if P n+ i(zo) = 0, then from (|3.4|) it follows that zo is an 
eigenvalue of J and the vector eo = (1, P 2 (zo), . . . , P n (zo)) T is corre- 
sponding eigenvector. 

Let z be a root of P n+ i(z) of the multiplicity I. Then 

P n+ l(z ) = P' n+ l{*«) = ■ = Pn+l(^) = 0, ^|l(^0) ± 0. 

Differentiating (|3.4|) we get 
(3.8) 

' hP[(z) + ai P^z) = zP[{z) + P 1 (z) 

a x P[{z) + b 2 P' 2 {z) + a 2 P.^z) = zP^z) + P 2 {z) 
< 

a n - 2 P' n _ 2 (z) + 6n_li^_!(z) + CLn^P'^z) = zP' n _ x {z) + P n _l(z) 

^ a^P^z) + 6 n P^z) + a n P' n+l {z) = zP' n {z) + P n (z) 

Recall that P\(z) = 1 and a n = 1. Substituting z = z we obtain that 
the vector e x = (0 P 2 {zo) ■ ■ ■ P^(zo)) satisfies (J — z I)ei = e . 
Continuing the differentiation of ([3.8)1 . we prove that the vectors given 
by ()3.6|) satisfy the equalities (J — z -0 e fc+i = e^, for = 1, . . . , Z — 
2. This means that the number zq is the root of the characteristic 
polynomial det(J — zT) and the multiplicity of this root is greater or 
equal I. 

Suppose now that the number zq is a root of the characteristic poly- 
nomial det(J — zl) of the multiplicity I > 1. Then as was proved 
above the number zq is the simple eigenvalue of J, zq is a root of the 
polynomial P n+ i(z) and eo is the corresponding eigenvector. Suppose 
that x = (x±, x 2 , . . . , x n ) T satisfies the equation (J — ZqI)x = 6q. Since 
the vectors x + Ae satisfy the same equation for arbitrary A, we may 
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assume that x\ = 0. Then we obtain the linear system 

' ai x 2 = P1O0) 
b 2 x 2 + a 2 x 3 = zx 2 + P 2 (z ) 

Pi-l 

(*b) 

Comparing the last linear system with (J3.8|) when z = Zq and taking 
into account that P[{z) = for all z, we get that x 2 = P 2 (z ), . . . , x n = 
P' n (zo) and P^ +1 (z ) = 0. So, x — e±. Following the same way as above, 

we get that P" +1 (z ) = ■ ■ ■ = Pn+i\ z o) — 0, i.e. the multiplicity of z 
as a root of the polynomial P n+ i(z) is greater or equal /. 

Thus, we proved that the roots and their multiplicities of the poly- 
nomial P n+ i(z) coincide with the roots and their multiplicities of the 
characteristic polynomial det(J — zl). 

The same is true for detf^Jji^j — zl) and for the polynomial P/ c+ i(z), 
k — 0, 1, . . . , n. Taking into account ()3.5|) . we get that ()3.7|) holds. □ 

Proposition 3.2. Let J be an n x n complex Jacobi matrix ^ for 
all k. Then the matrices J and J[i, n -i] have no common eigenvalues. 

Proof. Suppose that z is an eigenvalue of the matrices J and J[i tn -l\- 
Then by Proposition EH] we have P n (z ) = P n+ i(z ) = 0. Since 

fln-i-Pn-i(^o) + b n P n (zo) + a n+ iP n+ i(zo) = ZoP n (zo), 

we get P n -i(z ) = and from ()3.4|) it follows that P n - 2 (zo) = ■ ■ ■ = 
Pi( z o) — 0. But P\(z) = 1 for all z. Contradiction. Thus, the matrices 
J and J\i t n-i] have no common eigenvalues. □ 

3.3. Semi-infinite Jacobi matrices. In the following we will con- 
sider bounded Jacobi semi-infinite matrices of the form ()3.2j) . For the 
corresponding Jacobi operator 

the number z G C is an eigenvalue of J •<=>- 

(3 ' 9) £|P fc (z)| 2 <oo 

k=l 

and the corresponding eigenvectors are A (P 1 (z), P 2 (z), . . .) T , A G C. It 
was established in [7j that 

/;(./) <!zGC:sup ELil-W| 2 

^ n>x \a n \*[\P n {z)\ 2 + \P n+l {z)\ 2 ] 
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The function 

(3.10) mj(z) = ((J-zI)- 1 5 1 ,5 l ) 

is called the Weyl function of J. Note that the definition of the Weyl 
function in the form (|3.1U|) for a non-self-adjoint Jacobi matrix is given 
in [6 . The Weyl function has the following Taylor expansion at infinity 

™M = -L, z i+i 

1=0 

and the continued fraction (J-fractions) expansion |Mj 



(3.11) m,j{z) 



T — a 1 —a 2 a n-i 



; - bx + z - b 2 + z - b 3 + . . . + z — b n + . . . 
It follows that 



1 h bf + ai „ ( 1 

mj(z) J^A+ol- 

z z l z A V ^ 4 



oo. 



3.4. The m-function approach in inverse spectral problems for 
self-adjoint Jacobi matrices. It is well-known ^Hj that a finite real 
n x n Jacobi matrix J has n distinct real eigenvalues. A bounded real 
semi-infinite Jacobi matrix determines a bounded self-adjoint Jacobi 
operator in ^(N) with the simple spectrum. If E(t) is the orthogonal 
resolution of identity for the operator J, then dp(t) = d(E(t)5i, Si) is a 
probability measure supported at n points in the case of an n x n Jacobi 
matrix and on a finite interval of the real axis in the semi-infinite case. 
The probability measure dp is called the spectral measure of J. The 
corresponding Weyl function 

f dp(t) 
mj(z) = )—■ 

belongs to Herglotz-Nevanlinna class. 

The inverse spectral problems for finite and semi-infinite self-adjoint 
Jacobi matrices were studied in [IH], |H|, EH, E2, E2- The fol- 
lowing theorem is established by M. Stone for possibly unbounded 
self-adjoint operator. 

Theorem 3.3. [I], Let A be a self-adjoint operator with simple 

spectrum in a separable Hilbert space TC. Then there exists an orthonor- 
mal basis in 7i in which the matrix of A is a Jacobi matrix with the 
conditions 



(3.12) 



bk are real numbers, at > for all k 
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The construction of the corresponding self-adjoint Jacobi matrix can 
be provided by orthogonal normalization of l,x, x 2 ,... with respect 
to the measure dp or by means of m-functions (see [15J). We briefly 
describe the later approach for finite matrix and will keep notations of 
[To] . Denote by J\k, n ]i k = 2, . . . ,n the Jacobi matrix obtained from J 
by deleting k — l top rows and k — 1 left columns of J, J[i, n ] — J- Let 5k 
be A; x 1 column 5k = (0 ... 1 ... 0) T . Define the functions 

(3.13) m + (z, k - 1) = ({J[k,n\ - ziy l 5 k , 5 k ) , k = 1, . . . , n - 1. 

Thus m + (z, 0) = mj(z). As in [!!>] we will use the notation m + (z) for 
the function mj(z). 

Since the functions m + (z, k — l) has the expansion in a neighborhood 
of infinity 



m+(z,k- 1) = - 



]5k, 5k 



z i+i 



from (j3.1j) it follows that 
(3-14) 

m 4z,k-l)~- 1 -- b ±- b l±^ + o(±), * = l,...,n-l, 

lb b 2 / 1 
m + (z,n-l)~-±-%-% + ^ 

and (|2.16|) yields the following relations 

(3.15) a 2 k m + (z,k) + — — = b k - z, k = l,...,n. 

m + (z, k — l) 

Theorem 3.4. Every n-point probability measure is the spectral mea- 
sure of a unique n x n Jacobi matrix. 

Proof. We will follow the approach considered in ^3]. Let dp be a 
probability measure supported at n real points tx,...,t n . Define the 
Herglotz-Nevanlinna function 

f dp(t) -A a k 



t — z tk — z 

k=l K 

n 

where ak > and ^ = 1. Then in the neighborhood of infinity one 

fc=i 

has 

(.2 i „2 



- - +0 ' 



2! Z 2 Z 3 V Z 4 
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where 

n n 

h=} J ac k tk, b\ + a\ = 2J OL k t\ . 

k=l k=l 

One can find also 61 and a x > by 



1 



b\ = — lim 2 2 ( m(z) H — 
2^00 y 2: 

/ \ b 2 
a\ = -b\ - lim z 3 I 771(2;) H h 

^^oo y z z 1 

Define the Herglotz-Nevanlinna function m + (z, 1) by 

. , 1 / 1 
m + (z, 1) = — I 61 - z 



m(z) J 



This function has an expansion 



m+(z,l) r 5 — + 



Z Z 2 Z 3 » ' 4 



in the neighborhood of infinity. Find 62 and a 2 > by means of limits 
at infinity and define m + (z, 2) similarly Continuing this process up to 

m + (z,n - 1) = ( 6„_i - z 



"n-l \ 
and finding 6 n by 

1 



m + (z,n — 2)J b n — z 



b n = — lim 2 rn + (z, n — 1 

^^oo y ' 2; 

we obtain the real numbers {bk} k= i and the positive numbers {dkjkZi- 
Let us construct a self-adjoint Jacobi matrix J of the form (j3.1|) . Then 
the m-function 

m+(*,0) = ((J-z/)- 1 ^,^) 
coincides with m(z). □ 

Remark 3.5. TTie following J -fraction expansion of the function m + (z) 

(3.16) m + (z) = ^- -=4" Z^i 

z-6i+2-6 2 + 2- 63 + ■ ■ ■ + z-b n 

holds. 



For semi-infinite Jacobi matrices with conditions ()3.12|) and (jl.3|) the 
approach in [T5J is based on the representation of m(z) as a continued 
J-fraction of the form (|3.11|) and the following result: 
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Theorem 3.6. 15 . Suppose that m(z) = f(t — z)~ 1 dp(t), where dp 
is a probability measure on [— C, C] whose support contains more than 
one point. Let 



h = J dp(t), a{ = I t z dp{t) - bf 

and let 

mi(z) = \ ( &i 



m(z) 

Then m\{z) = J (t — z)~ l dp\(t), where dpi is a probability measure also 
supported on [-C, C] . 

We give one more application of the m-function approach which we 
will use later. 

Proposition 3.7. Let m(z) be a Herglotz-Nevanlinna function which 
is odd in some neighborhood of infinity and lim 2 _, 00 zm(z) = — 1. Then 
the diagonal entries of the corresponding self-adjoint Jacobi matrix are 
equal zero. 

Proof. Since m(—z) = —m(z), the function m(z) has the following 
Taylor's expansion at infinity 



k=0 



( 3 - 17 ) ™(*) = -EiSTT'ft = 1 - 



Let J be a bounded self-adjoint Jacobi matrix whose Weyl's function 
coincides with m(z). Then due to (|3.17|) we obtain b\ = 0. By (|3.15|) 

we have 

m+(z,l) = -j ( -z 



m{z) 

It follows that m+(— z, 1) = —m+(z,l). Now ()3.14j) yields 62 = 0. 
Hence by induction and again using ()3.15j) and (|3.14jh we obtain that 
b k = for all k. □ 

3.5. Mixed given data and uniqueness for finite self-adjoint 
Jacobi matrixces. Let J be an n x n self-adjoint Jacobi matrix. Fol- 
lowing [T3] we will consider 6's and a's as a single sequence {cfc}^ 1 , 
where c^k-i = b k and c^k = a&- The next theorem is established by 
H. Hochstadt in |2Tj . 

Theorem 3.8. Let J be annxn self-adjoint Jacobi matrix. Sup- 
pose that c n+ i, . . . , c 2n _i are known as well as the eigenvalues z\, . . . , z n 
of J. Then ci, . . . , c n are uniquely determined. 
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F. Gesztesy and B. Simon in [T5| proved the following generalization 
of the Hochstadt theorem. 

Theorem 3.9. ^H]- Let J be an n x n self-adjoint Jacobi matrix. 
Suppose that Cj + %, . . . , c 2n _i are known as well as j of the eigenvalues. 
Then C\, . . . ,Cj are uniquely determined. 

Note that it is not necessary to know which of the j eigenvalues 
one has, and there may be no matrix consistent with the data as in 
Theorems EE and EH (see [TJ ) . 

Since later on we will consider a non-self-adjoint analog of the Hochstadt 
and Gesztesy-Simon theorems, let us mention the key moment of the 
proof of Theorem 13.91 in [To] . Apart from m + (z, k)- functions the fol- 
lowing m_ - functions are used in [15J: 



(3.18) m-(z, k) = ((J[i,&-i] - zl) 4_i), k = 2, 3, . . . , n. 



It is established in [15] that for any eigenvalue Xj of a self-adjoint nxn 
Jacobi matrix J and for any k — 1, . . . , n the equality 



holds, where the equality in (|3.19J) includes the case that both sides equal 
infinity. 



4. Dissipative Jacobi matrices with a rank one imaginary 

PART 

4.1. Finite and semi-infinite dissipative Jacobi matrices. In the 

following we will consider the Jacobi matrix of the form ()3.1|) and (|3.2j) 
but now we suppose that condition (jl.4|) is fulfilled. 

Proposition 4.1. Let J be a non-self-adjoint n x n or semi-infinite 
Jacobi matrix of the form f!3.1|) with conditions (jl.4j) and (jl.3|) . Then 
J defines a prime dissipative operator with a rank one imaginary part 
in C n (^(N)) and its Livsic characteristic function takes the form 



(3.19) 




(4.1) 



W(z) = 1 - 2ilm 6i (( J - zl)- x 8 u 5 X ) . 



Proof. Let J* be the adjoint matrix to J and let 



J R = \{J + J*), J i = ^ J ~ r) 
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be Hermitian components of J. One has for n x n case 
(4.2) 



J, 



/Re 61 


a 1 








• 


■ \ 




/Im b\ 








ai 


h 


a 2 





• 




















h 


o 3 


• 




, Ji = 





















a n-l 










V • 








■ a n _ x 


K I 




V 









and for semi-infinite case 
(4.3) 

/Re 61 



J, 







a 1 

«2 





63 






«3 









V 



7 



/imfoi 











V 



and for every x G 

where 
(4.4) 



: n (z 2 (N)) 

2J lX 



fag) g, 



g = a/2 Im 61 5i . 
By (12.5)1 one obtains (|4.1jl . From it follows that J is prime. 



□ 



Remark 4.2. One can easily show that a non- self- adjoint finite or 
semi-infinite bounded tri-diagonal matrix A with diagonal entries {dk}, 
sub-diagonal entries {ck}, and super- diagonal entries {c~k} has a rank 
one imaginary part if and only if there exists a diagonal nonzero entry 
d m such that 

17 I |2 

and remaining entries equal zero. It follows that the Jacobi matrix J 
with the conditions Im&i = • • • = Im6 m _i = 0, Im6 m > 0, Im6 m+ i = 
• • • = ; ak > for all k is also dissipative with a rank one imaginary 
part. But in this case such a matrix might not be prime. For example 
the 3x3 matrix of the form 

/ 1 ai ' 
J — ( a i —l + it CI2 
\0 a 2 1 

is dissipative for t > and 

ran J/ = {X5 2 , A G C}. 
But span{^2, JS 2 , J 2 S 2 } 7^ C 3 . So, this matrix is not prime. 

The next proposition easily follows from ()2.6|) . (j2.16j) . and ()3.10|) . 



A 
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Proposition 4.3. Let J be a bounded non-self-adjoint finite or semi- 
infinite Jacobi matrix of the form (J3.1|) with condition (jl.4j) . Then 
for the Weyl functions mj(z) and mj R (z) of J and Jr and for the 
characteristic function W(z) there are the following relations: 



™j(z) = (W(z) - 1) 



_ i W(z) - 1 
mjR[Z) ~ Imb, W{z) + 1 

mj(z) 



%m jR \z) 



m J R { z ) 



i — Im fox rrij R (z) 
imj(z) 



i + Imfox mj(z) 
Moreover, the function 



mj(z) - mj(£) + 2zImb 1 mj(z)m J (£) 
&j{z,£) := = 

is a nonnegative kernel, the function 

1 

mj(z) 

is a Herglotz-N evanlinna, and at infinity 



1 \ 1 fo 2 bj + ai fl 

+ h-z) r ^- + O — 



af \ mj(z) ) z z 2 z 3 \z 4 / 

T/ie points of a non-real spectrum of J are solutions of the equation 

i 

m J R \ z ) = t — T~i lmz > °- 
imox 

From PropositionOl Theorem \F7\ and formulas (EH, (ETTUJ) . (|2~T2l . 
(I2.1HJ) . and (jrS.lOj) we get the following statement. 

Proposition 4.4. Lei J be a non-self-adjoint nxn Jacobi matrix of the 
form (|3.1|) with conditions fjl.4|) and fe£ £ 1; z 2 , • • • , z n be the eigenvalues 
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of J counting algebraic multiplicity. Then the following formulas 

n 

z- Z k 



W(z) = H 



. z- z k 

k=l 

n n 

n o - z$ - no- z k ) 

I \ 1 k=l k=l 



( 4 -5) ' " 2Im61 f[(z-z k ) 

k=l 

n n 

Y[(z-zt)- EI ( z - Z k 

I \ % k=l k=l 



lmbl f[{z-zi;) + f\{z-z k ) 



k=l k=l 



hold. 



Proposition 4.5. Let J be an n x n Jacobi matrix, n > 3, with the 
conditions (jl. 4)1 . Then the matrices J and Jh, n -2] have no common 
eigenvalues and the relations 

(46) _f^^__^ Im p<^) (0) = Im(ll 

hold for all k — 1, . . . , n, where the polynomials {Pk{ z )}k=i are defined 

Proof. Because all J\im are dissipative and prime, their eigenvalues 
belong to C + . Suppose that z is an eigenvalue of J and J[i, n _2]- Then 
by Proposition !3.1I P^„i (zn) = P n+ i(z ) = 0. From ()3.4|) we have 

a n _iP n _i(z ) + {b n - z )P n (z ) = 0. 

Since Im^o > and b n is real, we get that P n (zo) = 0, i.e. z is a 
common eigenvalue for the matrices J[i >n _2] an d Jn n _i]. By Proposition 
13.21 this is impossible. 

Let {zj}j =1 be the roots of the polynomial P k+ i(z). Then Imzj > 0. 
From ()3.7|) it follows that 

■n 

J2 z i = spj[i, fc ]. 

But ImSpJ[i 5 fc] = Imbi. On the other hand 



^ Q 1 Q 2 • • • flfc p(fe-l)/ n x 



Therefore, (JOf holds. □ 
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If J is a semi-infinite bounded Jacobi matrix of the form (|3.2|) and 
with conditions (jl.4|) . then by Proposition 14.11 it determines a prime 
dissipative operator with rank one imaginary part. Therefore its non- 
real spectrum is a denumerable set of eigenvalues of finite algebraic 
multiplicities and from ()3.9|) 

*p(J) = |^eC:^|P fc (^)| 2 <oo|, 

where {Pk(z)} are corresponding polynomials defined by ()3.4|) . By 
Theorem 12.31 the closure of the linear span of all corresponding root 
subspaces coincides with ^(N) if and only if 



Im - 



Im&i 



(4.7) 



J, 



Example 4.6. Let J\ be the semi-infinite Jacobi matrix of the form 

/ il 1/2 0- 
1/2 1/2 0- 
1/2 1/2 0- 

V 

where / > 0. The real part is the Jacobi matrix 

/ 1/2 0- 
1/2 1/2 0- 
1/2 1/2 0- 



7 



(4.8) 



\ 



7 



It is known [S] that the Weyl function of Hq is the function 



m Ho (z) 



2(Vz^l-z), Z eC\[-l,l]. 



The spectrum of Ho is continuous and coincides with the interval 



[—1,1]. Therefore, ||f/o 
tion 



1. The Weyl function of Ji is the func- 



mjXz) 



l-z) 



The non-real spectrum of J; we find solving the equation 



Vz 1 



i 

21' 



This equation has a unique solution from the open upper half-plane 

.4/ 2 -l 



4 
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for I > 1/2 and has no solutions for I G (0, 1/2]. Thus if I E (0, 1/2] 
the spectrum of J coincides with [—1,1] and if I > 1/2 the spectrum 



of Ji is {i(Al 2 
Note that 



1)/4}U[-1,1]. 
1 



lim 



lim 



sf-i m Ho (z) xii m Ho (z) 2 
Results from [23], jl], [3] yield that the Jacobi matrices of the form 

/ x 1/2 ■ ■ -\ 

1/2 1/2 
1/2 1/2 

V 

determine contractive Jacobi operators in ^(N) if and only if \x\ < 1/2. 

Example 4.7. Let J\ be the semi-infinite Jacobi matrix of the form 

/ il l/y/2 • • A 
l/y/2 1/2 0- 
1/2 1/2 0- 



• 7 



(4.9) 



Ji 







V 



7 



where / > 0. The real part of J; is the Jacobi matrix 

/ l/y/2 • • -\ 
l/y/2 1/2 • • 
1/2 1/2 • • 

v / 



(4.10) 



Hn 



It is known [H] that the Weyl function of Ho is the function 

1 



m 



Ho 



» 



y/7< 



«eC\[-i,l], 



The spectrum of i?o is continuous and coincides with the interval 
[—1,1]. Therefore \\Hq\ \ = 1. The Weyl function of is the function 



mj/z) 



iy/z^ 



l + l 



The non-real spectrum of J\ we find solving the equation 

y/z 2 - 1 = il. 

This equation has a unique solution from the open upper half-plane 
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for I > 1 and has no solutions for I G (0,1]. Thus if I G (0,1] the 
spectrum of Ji coincides with [—1, 1] and if I > 1 the spectrum of Ji is 

0VF^i}u[-i,i]. 



Note that 



lim |m& (x)\ =lim|mfr (x)\ 

xT-l H ° x|l 1 H \ J\ 



+00. 



From |22] it follows that the Jacobi matrix Hq is a unique contractive 
matrix in ^(N) among Jacobi matrices of the form 



( x l/y/2 • • \ 
1/V2 1/2 • • 
1/2 1/2 • • 

V J 



4.2. The m-functions of a finite dissipative Jacobi matrix with 
a rank one imaginary part. The next proposition establishes the 
analogues of the relations (J3.19|) for a dissipative Jacobi matrix. 

Theorem 4.8. Let J be an n x n Jacobi matrix with conditions (11.4)1 . 
Letmj(z) be the Weyl function of J defined by (|3.10|) and let the func- 
tions m + (z, j) and m_(z, j) be defined by (|3.13|) and (|3.18|) . If Zq is the 
eigenvalue of J of the algebraic multiplicity I then for all j = 1, . . . , n 
the following relations hold 



(4.n; 



i 



m-{z,j + 1, 



(p) 



0. 



z=z 



Proof. Represent the matrix J in the block-matrix form 

J -- 



A 2 i 



A 12 

J[j+l,n] 



where A 12 is j x (n — j) matrix of the form 



A 



12 



/o 




a ; 











o \ 








and An. = A* 12 (is (n-j) x j matrix). Let Hi = C j . Then from (I2~TBT) 

x r^i = (J M -zl- A l2 (J [j+hn] - zi)- i A 2l y l 



it follows that 
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An{J\j+i, n ] ~ zl) 1 Mi 



Now from ()2.8j) we get 
(4.12) 



((J-zI)-%5 J ) = T 



( 


••• 








••• 








••• 





v° 


••• 







>J+1) 










a 2 m+(z,j)m^(z,j + 1) 



a j m + (z,j)J 



j = 1,2,..., 71. 



Let j > 2. By Cramer's rule we have 

det(J [ i ij _i] - zl) det(J b+ i in] - zl) 



and 



det(J- zl) 
det(J[ij_i] - zl) 



det(J[ij] — zl) 
for j = 2, . . . , n. From ()4.12j) it follows that 
det(J [i+ i in] - zl) _ 1 



det(J — zl) det(J[ij] — zl) — a 2 m + (z,j) det(Jjij_i] — zl) 

Let Zq be an eigenvalue of J of the algebraic multiplicity /. Because J is 
dissipative, we have Imz > 0. Note that J\j+i, n ] is self-adjoint Jacobi 
matrix, the Herglotz-Nevanlinna function m + (z,j) is holomorphic at 
zq and lmm + (zo, j) > 0. It follows that the function 

det(J [i+1)TO ] - zl) 

det(J - zl) 

has the pole at z of order I. Therefore 

det(J [i+M -zl) 
hm (z — Zq) p TT77 — — T\ — = oo, p = 0, 1, . . . I — 1 



and 



det(J-zI) 
lim(s-*b)' dBt ^ + !' n] ~ Zl) is finite. 



2^2 



det(J- zl) 

(z - z ) p 



This yields 

(4.13) lim - — ^ — — — r = oo 

2^20 det(J[ij] — zl) — aj7n + (z,j) det(Jjij_i] — zl) 

for p = 0, 1, . . . , I — 1. In particular, for p = we obtain 

(4.14) det(J[ij] - z I) - a 2 m + (z ,j) det(J[i d _i] - z I) = 0. 
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By Proposition 13.21 the matrices Jmj-ii and J\ij] have no common 
eigenvalues and since m + (z ,j) ^ 0, (|4.14j) implies that 

det(J[ij] - z I) ^ 0, det(J[ij_!] - z I) ^ 

and 

1 det( Jh ,-i — zqI) 

4.15 , — - = \ M °> = a 2 jm+ (z ,j). 

m_(z ,J + l) det(J[ij_i] - zqI) 3 

In some neighborhood of z we can rewrite ()4.13|) as 

lim — k t r = OO 

*-«o det^u-i] - zl) {m-(z,j + l)) 1 - a^m + {z,j) 

for p = 0,1, — 1. It follows that the number zq is a zero of the 
function 

a 2 j m + {z,j) 



m_(z,j + 1) 3 

of order I. Thus, we get (14.11)1 for j > 2. 

Let j = 1. Since m_(z, 2) = (b\ — z)^ 1 , relation (|4.12j) takes the form 

m j( z ) = 7 \ — r^V 

b\ — z — aim + {z } 1) 

If Zq is the eigenvalue of J of the algebraic multiplicity I, then due to 
(14. 5|) the Weyl function mj(z) has a pole of order I at zo- I n the same 
manner as the above we get that zq is a zero of order I for the function 

h-z- afm + (z, 1) = ^— - a\m + (z, 1). 

Therefore, (|4.11|) holds and for j = 1. □ 

Since the integer j is an arbitrary, we obtain as a by-product the 
following corollary. 

Corollary 4.9. Let J be an n x n Jacobi matrix with the conditions 
(|1.4J) . Then the matrices J\ij\ and J have no common eigenvalues for 
every j = 1, 2, . . . , n — 1. 

5. Inverse spectral problems for finite Jacobi matrices 
with a rank one imaginary part 

5.1. Reconstruction of a finite dissipative Jacobi matrix from 
its eigenvalues. Let J be a non-self-adjoint n x n Jacobi matrix of 
the form (|3.1|) satisfying conditions (jl.4|) . Then by Proposition 14. II the 
corresponding operator in C™ is a prime dissipative operator with a 
rank one imaginary part. Therefore the matrix J has only non-real 
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eigenvalues with positive imaginary parts. The next theorem estab- 
lishes that an arbitrary n non-real numbers counting multiplicity from 
the open upper half-plane determine uniquely some dissipative n x n 
Jacobi matrix with a rank one imaginary part. 

Theorem 5.1. Suppose that z\, . . . , z n are not necessarily distinct com- 
plex numbers with positive imaginary parts. Then there exists a unique 
nxn Jacobi matrix with entries satisfying conditions (jl.4|) whose eigen- 
values counting algebraic multiplicity coincide with {zk}k=i- 



Proof. Let 

(5.1) 

Then 

Let 

Define 
(5.2) 



w(z) = n 



z k 



k=l 



lim z(W(z) 



2i^Im 



Zk- 



k=l 



k=l 



m + {z,0) 



i W(z) 



c W{z) + 1 

The Herglotz-Nevanlinna function m + (z, 0) has the expansion at infin- 
ity 



m+(z,0) 



1 

z 



b b 2 + a\ 



+ 



and determines a probability measure supported at some n points. By 
Theorem 13.41 there exists a unique self-adjoint nxn Jacobi matrix 



H 



/ b a\ 
di b 2 
a 2 





a 2 
b 3 






a 3 



\ 



V 



a n _i b n ) 



satisfying conditions (j3.12|) such that 

m + (z,0) = ((H-zI)- 1 5 1 ,5 1 ). 
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J = 
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fb + ic Oi 

CL\ 62 O2 

a 2 63 «3 



V 



0"n-l 

a„_i 6 n / 



/c 





V 



•0 0, 



It follows that 

2J r x={x,g), x e C™, 

where g = \[2cb\ £ C n . 

By Proposition 14. II for the characteristic function of J we obtain 

w(z) = 1 - 2i c ((J - z/)- 1 ^, 8x) . 

Its fractional-linear transformation 

. w(z) — 1 



VIZ 



w(z) + 1 



has the form 



v{z) = - ({H - zl)~ l g, g) = c ({H - zI)~H x , 5 ± ) = cm + (z, 0). 



From (|5.2j) we get 



viz) 



W{z) 



1 



W(z) + 1 

Therefore w(z) = W{z). By Theorem 12 . 71 the eigenvalues of J counting 
algebraic multiplicity coincide with {zf.}. □ 

Example 5.2. Let us construct 3x3 dissipative Jacobi matrix with 
eigenvalues Zi — % of multiplicity 2 and 2; 2 = 2i (of multiplicity 1). By 
Theorem 12.61 the characteristic function of J is 

/ y — 1 \ 

W(z) 



Then 



z + i\ z + 2i 



z — i J z — 2i 



lim z(W(z) - 1) = 4z. 

2— >oo 

Note that 21m zi + Im 02 = 4. Let 



m + (z,0) = - 



1 _ i (z + i) 2 (z + 2i) - (z - if(z - 2i) 
4 W(z) + 1 ~ 4 (z + i) 2 (z + 2i) + (z - i) 2 (z - 2i) 
-2z 2 + 1 
2z 3 - lOz ' 
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Then 



m(z) 



z — 



9/2 



z — 



1/2 



It follows that bi = 4i, b 2 = b 3 
Jacobi matrix takes the form 



J 



0, ai = y/9/2, a 2 = y/T/2. The 




Example 5.3. In order to construct the dissipative nxn Jacobi matrix 
with the eigenvalue z = x + iy of algebraic multiplicity n (y > 0) it 
is sufficient to construct the nxn Jacobi matrix J n with the eigenvalue 
z Q = i of algebraic multiplicity n. The characteristic function of such 
a matrix is 

w n (z) = ( £±i' 

\z — I 

The corresponding Weyl function of the real part is 

% (z + i) n - (z - i) n 



m n (z) 



n (z + i) n + (z — i) 



Because SpJ = in, we get Imfei = n. Expanding the function m n (z) 
via continued fraction we obtain the real part of J n . In particular 



m 2 (z) 



-1 



z — 



m 3 (z) 



-1 



z — 



8/3 



1/3 



m±{z) = 



z + 



-1 



-5 



z + 



-4/5 



z + 



-1/5 



m 5 (z) = 



z + 



-1 



-8 



z + 



-7/5 



z + 



- 16/35 



z + 



Therefore 



J 2 



2i 
1 



J 3 



/ 3i 


2\/2 
V3 


o\ 







1 




V3 


V 


1 

V3 
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J 4 



( Ai 




V5 


2 







2 



1 







f 5i ^ 



i 




/ 



V 



2^ 







2V2 
VE 



VI 

Vr> 








VI 
Vf> 



4 








4 


1 








1 





/ 



Note that the eigenvalues of the real part of J n are solutions of the 
equations 



z + i 



that are real numbers 

X k = — cot 



71 + 2nk 
2n 



-1, 



k = 0,1 



n — 1. 



5.2. The connection between triangular and Jacobi forms of 
a prime dissipative operator with a rank one imaginary part. 

Let Zi, Z2, ■ ■ ■ , z n be not necessarily distinct complex number with pos- 
itive imaginary parts. According to Theorem 15 . II there exists a unique 
n x n Jacobi matrix 





fh 


a 1 





• 






b 2 


tt? 


• 


J = 





«2 


63 


a 3 • 



V ■ ■ • -0 a n _i b n J 

with entries satisfying conditions ()1.4|) whose eigenvalues coincide with 
{zk\1=i counting algebraic multiplicity. On the other hand such a 
matrix by Theorem 12.71 is unitary equivalent to the triangular matrix 
of the form 

i/3i(3 n \ 

Z n ) 

where 

f3 2 

(5.4) z k = a k + i^ = lmz k , (3 k > 0, k = 1, . . . , n. 



(5.3) 



J 



fzi 




Z2 







NON-SELF-ADJOINT JACOBI MATRICES 
From (|5.Hj) it follows that 

/ Pi P1P2 ■ ■ ■ PiPn\ 
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21m J = 

and 

(5.5) 

where 

(5.6) 

Let U be unitary matrix 
U = 

such that 
(5.7) 



P1P2 Pi 



PA 



\PlPn fcPn ■ ■ ■ PI J 

21m J = (•,#) if, 



//3i\ 

ft 



\PJ 



{ Uu U12 



Ui n \ 
U2n 



UJ=JU 
Ug = g 



where g = \ / 2lmb\ 5\. Next we present an algorithm which allows to 
find the Jacobi matrix J and the unitary matrix U satisfying (|5.7|) . 
From (J5.7j) it follows 



(5.* 



UJ k g = (J) k g, k = l,2,...,n. 



Since U is the unitary matrix, we have ||g|| 2 
account (|4.4|) and ()5.6|) we get 



\g\\ 2 , and taking into 



\9\ 



21m b 1 ^Pl^Pl + ... + p 2 n = \\g\\ 2 . 



Thus 
(5.9) 



Im b\ = Im 



Zk- 



k=l 
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From Ug — g it follows 



(5.10) 



Vim z k 

Uki = —f= k = 1,2, ...,n. 



E Im Z j 



Relations (jSHJ) and (Ql yield 

Jg,g\ = (UJg,Ug) = (Jg,g). 
Taking into account (jl.lj) . (|4.4j) . (|5.3jl and ()5.6|) . we get 

/ 



/6iV 2Im 6i \ 
ax\/2 Im &i 




(5.11) 



V 



/ 



fc=2 



\ 



z2 + *E$)ft 

fc=3 



Therefore 



Since (Jg,g) 



(Jg,g) = 26ilm6!, 

n / n \ 

Jg, g ) , from (|5.4j) and ()5.9|) we get 



E Uj + * E PI] lmz i E lmz k 

= ^ ; — Re6i = 



E i m2; fc 
fe=i 



E i m ^ 

k=i 



Because U is unitary and UJg = Jg, we get \\Jg\ 
equality and (j5.11|) yield 



Jg 



|&i| 2 (2Im&i) +a?(2Im&i) = E 



E a 
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and hence 
/ 



«i 



n 




n 


2 


n 


E 


Zj + i 


E ft 


Im Zj 


E 




k 


=j+l 




3=1 



k=j+l 



2\ 1/2 



E Im2: fc 



V 



fc=l 



E 

fe=i 



Recall that by (|3.3|) 



(J m 1 g) m = a m _ia m _ 2 • • • ai\/2 Im&i 



and 



m+2 



(J m "^)n = 



for m = l,...,n. Let m > 2. Suppose that (J m 1 g)k are already 
known, where k = 1, 2..., m — 1. Then from we obtain 

/ 61(^-^)1 + a 1 (J w - 1 p) 2 \ 

ai(J m -^)i + hiJ^g^ + a 2 (J m ~ ^) 3 



a m _ 1 (J m -^) m _ 1 + 6 m (J m - 1 ( 7 )„ 

m 





It follows that 



m—l 



{J m g, J m ~ l g) Cn = J> m <7)i (J™- 1 *?); = £(J m <^ 
i=i j=i 



+ (a m _ 1 (J" l - 1 ^) OT _ 1 + 6 m (J m -^) m ) (J m "^) m . 
From 
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and (J m 1 g) m = a m~i a m-2 4 • • Oi\/2 Im&x we get the linear equation 
with respect to b rn \ 

m—l 



7m- 1. 



This equation can be solved since the coefficient (J m 1 g) m is nonzero 
by (|3.3|) . In order to find a m we use the relation 

2 



1^1 



which takes the form 

m—l 



/.7 1 



+ \a 

m—l 



(J) m g 



m—l 



+ 



Solving for a 2 m we can find a m . 

According to (j5.10J) the elements {uki}^ =1 are expressed by means of 
Zi, z 2 , ■ ■ ■ , z n . The equality UJg = Jg gives the following linear system 
with respect to the second column {u fc2 }fc = i of the matrix U: 

' u\\(Jg)\ + u 12 (Jg) 2 = (Jg)i 
u 2 i(Jg)i + u 22 (Jg)2 = {Jg)2 



, u n i(Jg)i + u n2 (Jg) 2 = (Jg)n 



Since {Jg) 2 = aw/I Im&i ^ 0, one can find {iik2}k=i- By induction the 
equality UJ m ~ 1 g = {J) m ~ 1 g enables us to find {uk m }k=i f° r fn < n. 

Finally the relation J = U~ 1 JU and already known entries {ak}2Z\ 
and {bk}2Zi allow us to find the entry b n . 

In conclusion we give formulas for the reconstruction of 2 x 2 dissi- 
pative Jacobi matrix J with a rank one imaginary part from its eigen- 
values zi, z 2 , and the corresponding unitary matrix U: 



J 





Re z\— Rc Z2 
Im zi+Im Z2 



1 I Imz! Im z 2 



Re z\ Im .Z2+RC Z2 Im z\ 
Im zi+lm Z2 
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u 



Im z\ 



Im Z\ + Im Z2 



Im z 2 



Im Z\ + Im Z2 



Im Z2 



Im 0i 



\ V Im Z\ + Im ^2 V I m z i + Im ^2 / 

Remark 5.4. T/ie 'presented algorithm of reconstruction of the unique 
dissipative non- self- adjoint Jacobi matrix with a rank one imaginary 
part having given non-real numbers as its eigenvalues allows to recover 
the set of tri-diagonal matrices with the same non-real eigenvalues. 

5.3. A non-self-adjoint analog of the Hochstadt and Gesztesy— 
Simon theorems. The next theorem is a non-self-adjoint analog the 
of Hochstadt |2T] and Gesztesy -Simon [TH] results (see Theorems 13.81 
andlHH)- 

Theorem 5.5. Let J be an n x n Jacobi matrix with conditions (jl.4|) . 
Suppose that the eigenvalues z±, . . . , Zf. £ C + counting algebraic multi- 
plicity li, . . . , l k are known as well as 

b±, ai, 62, . . . , a„_ r _i, 6 n _ r , 

where 

r = h + k H h Ik- 
Then a n - r , b n - r+ i, . . . , a n _i, 6 n are uniquely determined. 

Proof. The sub-matrix J[ l ri _ r ] is known. Therefore, the m_-function 
m_(z, n — r + 1) is known. The unknown Jacobi sub-matrix J[ n - r +i,i] is 
self-adjoint and the corresponding Weyl function m+(z, n — r) belongs 
to the Herglotz-Nevanlinna class. By Cramer's rule we have that 



L r -i{z) 
L r (z) 



where L r _i(z) and L r (z) are polynomials with real coefficients of de- 
grees r — 1 and r, respectively. Let 



$(z) 

By Theorem 14.81 we have 



, m + (2;, n — r). 



m_(z, n — r + 1) 



(p) 



k, p — 0,...L 



Therefore the values 



j = p = 0,...L 



are known. Since = $(z) and Imz,, > 0, the values 



k, p 



0,...l r 



44 



YURY ARLINSKII AND EDUARD TSEKANOVSKII 



are also known. Let $>(z) be a ratio of two polynomials with real 
coefficients of degrees r — 1 and r such that 

$<p\ Zj ) = = $ (p) (^)j = l,...,k,p = 0,...l j . 

Then for the function ^(z) = $(z) — &(z) we get 

y<p\ Zj ) = *W(^) = o, j = 1, . . . , fc, P = 0, . . .l r 

It follows that 



k 

*(z)=^(z)l[( 

3=1 



where the function if)(z) has no zeroes at the points {zj}^ =1 and {z]}^ =1 . 
The function *&(z) is the ratio of two polynomials. The numerator is 
of degree less or equal 2r — 1. Since r = 1% + . . . + we get ^(2) = 0. 
So, $(2) is uniquely determined. Since 

ai- r = - lim z$(z), 

z—>oo 

we get the function m + (z,n — r) which uniquely determines the self- 
adjoint Jacobi matrix J[ n _ r+ i jn ], i.e. the entries 6 n _ r+ i, . . . , a n _i, b n are 
uniquely determined. □ 

5.4. Extension and refinement of the Hochstadt and Gesztesy— 
Simon uniqueness theorems for tri-diagonal matrices. Let J be 

the n x n tri-diagonal matrix of the form (jl.ljl with the conditions 



(5.12) 



Im b\ > 0, bk = bk, k = 2, . . . , n, 
dk > 0, k 7^ p, a p = 0. 



Recall (see Subsection I3.5|) that we consider a's and fe's as a single 
sequence {c^}^ , where c 2 fe_i = 6^ and c 2 fc = a^. The next theorem 
is the extension and refinement of the Hochstadt theorem [2Tj . 

Theorem 5.6. Let J be an n x n tri-diagonal matrix of the form (jl.lj) 
wii/i conditions (|5.12j) . Suppose that p non-real eigenvalues z%, . . . ,z p 
counting algebraic multiplicity and distinct real eigenvalues 
Zp+i, . . . , z n of J are given as well as c p+n+ i, . . . , c 2n _i. T/ien these data 
uniquely determine entries ci, . . . , c p+n . 

Proof. Since a p = the matrix J takes the block form 
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where Jn and J 22 are p x p and 
of the form 



n - 



p) x {n — p) tri-diagonal matrices 



(5.14) 



Jn 



CLi 




b 2 
«2 





«2 






o 3 



a p _i 



a p _i 6 P / 



(5.15) 



•I 



22 



'p+i 



bp+2 





Qp+2 



1 



V 



a n _i 



The matrix J22 is self-adjoint and therefore it has only real distinct n—p 
eigenvalues which are also eigenvalues of J. The matrix Jn is prime dis- 
sipative Jacobi matrix with a rank one imaginary part. It follows that 
Jn as well as J has p non-real eigenvalues counting multiplicity. By 
Theorem 15. li the entries ci, c 2 , . . . , c 2p _i can be reconstructed by means 
of the spectral data z\, 



z p . The distinct real numbers z p , 



Zri are 



eigenvalues of self-adjoint Jacobi matrix J22. By the Hochstadt theorem 
the entries c 2p +i, ■ ■ ■ , C2 P +( n -p) are determined uniquely by z p+ i, . . . , z n 
and by the known entries c p+n+ i, . . . , c 2n -i- D 

Next is the extension and refinement of the Gesztesy-Simon theorem 

EST 

Theorem 5.7. Let J be an nxn tri-diagonal matrix of the form (jl.lj) 
with conditions (|5.12j) . Suppose that p non-real eigenvalues Zi, . . . ,z p 
counting algebraic multiplicity and distinct real j eigenvalues of J are 
given as well as entries 

C 2 p+j+l, • • • , c 2n -i- 
Then these data uniquely determine the matrix J. 

Proof. From the proof of Theorem 15.61 and conditions (|5.12j) the tri- 
diagonal matrix J takes the form (J5.13j) . By Theorem 15. II the entries 
Ci, c 2 , . . . , c 2p _i can be reconstructed uniquely by means of the non-real 
spectral data zi,...z p . Using the Gesztesy-Simon uniqueness Theo- 
rem 13.91 the j distinct real eigenvalues and entries 02^4.^+1, . . . , C2 n -i 
uniquely determine the matrix J 22 . Therefore the given data uniquely 
determines matrix J. 

□ 
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Remark 5.8. Consider more general tri-diagonal n x n dissipative 
matrix J with ran Jj = span {5i} of the form 





fh 


«i 





• 






b 2 




• 


J = 





a-2 


&3 


a 3 • 



V • • • -0 a n _i b n / 

with the condition 

Im hi > 0, 6fe = bk, k — 2, . . . , n 
For such matrices the following statements are equivalent: 

(1) the matrix J has p non-real eigenvalues counting their algebraic 
multiplicity and p < n; 

(2) a p = and 7^ for k — 1, . . . ,p — 1. 

Actually, suppose that a p = and all % ^ for A; = 1, 2, ... ,p — 1. 
Then J takes the block form ()5.13|) . where J\\ and J 22 are p x p and 
(ra — p) x {n — p) tri-diagonal matrices given by (j5.14|) and ()5.15|) . The 
matrix J n is prime dissipative Jacobi matrix with a rank one imaginary 
part and J22 is self-adjoint. It follows that J has p non-real eigenvalues 
counting algebraic multiplicity. 

Conversely, suppose that the matrix J has p non-real eigenvalues 
counting algebraic multiplicity and p < n. Let J be the corresponding 
operator in C n . Since ran Jj = span {Si}, the subspace 

Hi = span {5i, JSi, . . . J n Si} 

reduces J, Jn = J\ Hi is a prime dissipative operator with a rank one 
imaginary part, J22 = J\ (C n QHi) is a self-adjoint operator. Hence Jn 
has only non-real eigenvalues that coincide with non-real eigenvalues 
of J. It follows that dim Hi = p. From ()3.3|) it follows that 7^ for 
k — 1, . . . ,p — 1, and a p = 0. 

6. A NON-SELF-ADJOINT ANALOG OF THE STONE THEOREM 

The next theorem is a non-self-adjoint analog of M. Stone's result 
(see Theorem 13. 3J) . 

Theorem 6.1. Let H be separable Hilbert space and let A be a bounded 
prime dissipative operator with a rank one imaginary part acting in H. 
Then there exists an orthonormal basis in H in which the matrix of 
the operator A is a bounded dissipative Jacobi matrix with conditions 

(Q. 
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Proof. Let dim 7-^ = n. Because operator A is prime, it has only n 
non-real eigenvalues {z k } C C + counting their algebraic multiplicity. 
By Theorem 12 . 71 the characteristic function of A takes the form 



w(z) = n 



k=l 



Z - Z k 



By Theorem 15.11 there exists a Jacobi matrix 

fbi ai 

ax b 2 a 2 
a 2 b 3 a 3 



V 



a„-i b n ) 



with conditions ()1.4|) whose eigenvalues coincide with {z^} and the 
corresponding characteristic function coincides with W(z). It follows 
that A is unitarily equivalent to the operator in C n determined by the 
matrix J. 

Let dimH = oo and let g e H such that 2Ajh = (h,g)g, h G H. 
Define the characteristic function of A 



W{z) = l-i({A-zI)- l g,g), z e p(A) 



It follows that 



Let c = | \g\ | 2 . Define 



lim z{W{z) - 1) 



i\\g\ 



V(z) 



W(z) 



and let 



Then 



W(z] 
m(z) 

m(z) 



\ = \{IAr-*I)- X 9>9) 



V(z). 



dp(t) 
t-z 



where the spectral measure has bounded support. Then there exists a 
unique Jacobi matrix 

/ b oi • • \ 

ai b 2 a 2 • • 
a 2 6 3 ffl 3 ■ ■ 



H 



V 



7 
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with real entries b, 62, b%, . . . and positive entries a±, . . . such that 
m(z) = ((H — zl)~ 1 5i, 5i) . Note that (jl.3|) holds because H defines a 
bounded operator in ^(N). Moreover, the entries of H can be found by 
means of continued fraction expansion 



m(z) 



"n-l 



Let 



b 3 + ... + z -b n + 



J 



\\9\ 
a 1 





2 /2 



a 1 
a 2 





a 2 
h 



/IMI72 



H + i 











03 

















V 



7 



The characteristic function of J is 

l-i|b|| 2 ((J-^)^i,5i). 

Since Jr = H, and m(z) = ((H — z/) -1 ^, 5i) , we get that the charac- 
teristic function of J coincides with W(z). Because matrix J is prime, 
the operator A is unitarily equivalent to J. Note that the entries of J 
can be also found using the continued fraction expansion 



M{z) 
where M(z) 



"n-l 



b x + z -b 2 + z - 63 



UW{z) - 1), and (3 



lim (iz( 

z— >oo 



+ z - b n + 
l-W(z))). 



□ 



From Theorem 16.11 it follows that the Jacobi matrices with condi- 
tions (jl.3|) and (jl.4j) provide new models for bounded prime dissipative 
operators with rank one imaginary part alternative to Livsic triangular 
models described in Subsection 12. < 



Remark 6.2. From Theorems \2.10\ and \6.1\ it follows that the equiva- 
lent statements (Hi), (iv), and (v) are also equivalent to the statement: 
The function M(z) is the Weyl function of some bounded Jacobi matrix 
with conditions (jl.4|) . 

Remark 6.3. In the recent paper ^3] it is proved that every cyclic self- 
adjoint operator in a Pontryagin space is unitary equivalent to some 
generalized Jacobi matrix 



NON-SELF-ADJOINT JACOBI MATRICES 



49 



The next theorem gives more information about the diagonal entries 
of the model Jacobi matrix. 

Theorem 6.4. Let A be a bounded prime dissipative operator with a 
rank one imaginary part acting in a separable Hilbert space 7i. Suppose 
that the characteristic function W(z) of A possesses the property 

(6.1) W(-z) = ^ y \z\>\\A\\. 

Then the corresponding to A (by Theorem \2. 71 ) bounded Jacobi matrix 
of the form (jrS.l)) with the conditions (jl.4j) possesses the property 

Re b x = b 2 = • • • = 0. 

Proof. Since (see (EHTA ) 

V(z)=i ™W- 1 l ,zep(A)n P (A R ), 

we have that in a neighborhood of infinity 

V(-z) = -V{z). 

Because 

where / is the positive eigenvalue of Ai, we get that the Weyl function 
of Ar is odd. Since mj R (z) = uia r (z), from Proposition 13.71 it follows 
that Re &i = = • • • = 0. □ 

Corollary 6.5. Let A be a bounded prime dissipative operator with a 
rank one imaginary part acting in a separable Hilbert space TC. Suppose 
that eigenvalues of A (counting algebraic multiplicity) are symmetric 
with respect to the imaginary axis and A has a complete system of root 
subspaces. If J is the corresponding to A bounded Jacobi matrix with 
conditions (|1.4j) . then J possesses the property Re b\ = b 2 = ■ ■ ■ = 0. 

Proof. By Livsic Theorem 12.61 the characteristic function W{z) takes 
the form 

N / - \ h N / i \h M i 

w(z) = TT ( i^i) TT I £±5t) TT LtSi, 

where {zj}f =1) {— Zj}f=i (J mz j > 0, Rezj < 0) are distinct eigenvalues 
of A, {lj}jL x are corresponding algebraic multiplicities, {iT]k}^=i, Vk > 
are not necessarily distinct pure imaginary eigenvalues of A, N and M 
are both finite when dim 7i < oo, and one of them or both are infinite 
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N 

when dim ft = oo, and M + 2 ]T \ j = dim H. Note that if M = 0, (i.e. 

there are no pure imaginary eigenvalues) then in the finite dimensional 
T~C its dimension is even. It is also possible that N = 0. 
It follows that 

w{ ~ z) = w(7r |2|>M ' 

By Theorem 16.41 we obtain that the corresponding Jacobi matrix pos- 
sesses the property Re b\ — b 2 = ■ • • = 0. □ 

Now we want to mention some connections of Theorems 15.11 16.11 and 
Corollary 16 . 51 with some result established by K. Veselic in [3H] and jHI] 
related to the spectral problems for the operator pencil of the form 

(6.2) z 2 M + zC + K, 

where M, C, and K are self-adjoint n x n matrices with real entries, 
M and K are positive definite, C is positive semi-definite. In [30 this 
problem is reduced to the eigenvalue problem for some 2n x 2n matrix 
A constructed by means of M, C, and K. The inverse spectral problem 
for ()6.2|) is considered in |3T] for the case 

(6.3) M = diag(mi,m2, . . . ,m n ), 



(6.4) 



and 



(6.5) 



K 



( h -h o o 

-&i £4 + k 2 -k 2 



C 



/7 





k n — 1 

1 kn—1 ~l~ kn/ 



000/ 



where mi = 1, {rrij}, {kj}, and 7 are positive numbers. In such situa- 
tion the following result is established. 

Theorem 6.6. [31J. Prescribe 2n eigenvalues from the open left half- 
plane which are symmetric (together with their algebraic multiplicities) 
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with respect to the real axis. Then there exist unique matrices M , K 
and C of the form (|6.3|) . (|6.4|) . (|6.5|) . respectively, such that 

(i) nij > 0, kj > 0, mi = 1, 

(ii) prescribed eigenvalues coincide with the spectrum of the operator 
pencil ()6.2|) . 

Note that in jHOj and jHI] it is supposed that the total number of 
negative eigenvalues counting their algebraic multiplicities is even. 

As it was mentioned by K. Veselic (in private communication) the 
matrix A is unitarily similar to In x 2n tri-diagonal matrix A 1 of the 
form 



/ 



with «, 



A' 



v^7 



-7 
-cti 





a i 


-ft 





ft 







«2 



-a r 



\ 



a, 
0/ 



in 



j- 



J 



,n, Pj = ^/k~J 



n — 1. 



The matrix J' = —iA' is a dissipative tri-diagonal matrix with a rank 
one imaginary part. Clearly, the matrix J' is unitarily similar to 2nx2n 
Jacobi matrix 



J 



/«7 cti 
c*i ft • 
ft a 2 o 



a, 



V 



o « n 0/ 



It follows that the inverse spectral problem considered by K. Veselic is 
reduced to the reconstruction problem of a dissipative Jacobi matrix 
with a rank one imaginary part from its eigenvalues. The symmetry 
property for the eigenvalues of the pencil ()6.2j) means that the eigen- 
values of J counting algebraic multiplicity are symmetric with respect 
to the imaginary axis. According to Theorem 15. II there exists a unique 
2n x 2n Jacobi matrix with entries satisfying conditions ()1.4|) whose 
eigenvalues counting algebraic multiplicity are prescribed numbers from 
the open upper half-plane and symmetric with respect to imaginary 
axis. By Corollary 16 . 51 the diagonal entries of the reconstructed Jacobi 
matrix possess the property 



Re&i = b 2 = • • • = b 2n = 0. 
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In conclusion we note that the proof of Theorem 16.61 in 31 differs from 
our proof of Theorem 15.11 and does not use the Livsic characteristic 
function. 



7. DlSSIPATIVE VOLTERRA OPERATORS WITH RANK ONE 
IMAGINARY PARTS AND THE CORRESPONDING SEMI-INFINITE 

JACOBI MATRICES 

Recall that a compact operator is called Volterra operator if its spec- 
trum consists of one point (zero). 

Theorem 7.1. [27 i Let 7i be a separable Hilbert space. Let A be a 
prime dissipative Volterra operator with a rank one imaginary part and 
let 



(7.1) Aih 

for every h, where I > 0, \\e\\ — 
the integration operator 



= l(h, e)e 

1. Then A is unitarily equivalent to 



(Ff)(x) = 2i / f(t)dt 



in the Hilbert space ^[0, 1]- 

In connection with Theorem 16.11 we will find a Jacobi matrix corre- 
sponding to Volterra operator. 

Theorem 7.2. Let A be a prime dissipative Volterra operator with 
a rank one imaginary part and let Ajh = l(h, e)e for every h, where 
I > 0, ||e|| = 1. Then A is unitarily equivalent to the operator in ^(N) 
determined by the dissipative Jacobi matrix 
(7.2) 



/ il 

i 



i 



J 



/1-3 




/1-3 

I 

/3S 





I 









I 









V 



g I 

v /(2n-3)(2n-l) ^/2n-l)(2n+l) 



Proof. As is well known j^ZI, OH] that if (|7.1|) holds, then the charac- 
teristic function W(z) — 1 — 2il ((A — zl)~ 1 e, e) of A takes the form 



W(z) = exp 



f) 
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It follows that 

V(z) = l((A R -zI)- 1 e,e) 



W(z) - 1 



tan - 



W(z) + 1 

By Proposition 12.21 the vector e is a cyclic for A R , therefore m(z) 
((An — z/) -1 e, e) is the Weyl function of A R . We get 



m(z) 



1 (I 
-- tan - 
/ V z 



We consider for simplicity the case I — 1. In this case m(z) = — tan (-) . 
We have to find a self-adjoint semi- infinite Jacobi matrix whose the 
Weyl function is — tan (~) . This function has the following continued 
fraction expansion [24J 



tanx 



x 

I 



X 

y 



X 

y 



X 

y 



X 

y 



X 



.. -2n + l - 

Comparing with ()3.11|) and taking into account that x = 1/z we get 

1 



0, at 



n 



1,2, 



(2n- l)(2n + l) 

It follows that the Jacobi matrix with the Weyl function 
the following 





tan (i) 



is 



/ o 

i 



Let 



yi-3 




/1-3 


1 




l 






1 

/5~7 









A 







^/(2n-3)(2n-l) y/2n-l)(2n+l) 



1 



Jo 



/1-3 




/1-3 


1 

/3~5 





1 

/3~5 









1 

7^7 









V 



1 o 1 

v /(2n-3)(2n-l) ^/2n-l)(2n+l) 



7 



Then the characteristic function of Jo is exp(j-). So the operator 
J is unitary equivalent to A with / = 1. For arbitrary / the Jacobi 
matrix J = U takes the form (J7.2j) and its characteristic function is 
exp(f). □ 
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Remark 7.3. The real part Tr of the integration operator (Ff)(x) 
i 

2i J f(t)dt takes the form 



(F R f)(x)=i\ I f(t)dt - I f(t)dt 
Its eigenvalues are 



21 



kez. 



(2k + l)vr ; 

From Theorem \7.^ it follows that these numbers are eigenvalues of the 
matrix 
(7.3) 



H 



m 
o 



A-3 





/3-5 







/3-5 





/5-7 



V 



i o ' 

• v /(2n-3)(2n-l) y , 2n-l)(2n+l) 



7 



Remark 7.4. zs known pQ, |2H| #iat t/jere are relations between the 
moments {7^ = J t k dp(t)}, 70 = 1 of a measure dp(t) and the entries 
{bk}, {ak\ of the corresponding self-adjoint Jacobi matrix 

(b x a 1 • • \ 
a\ b 2 a 2 • ■ 
a 2 63 a 3 • ■ 



J 



\ 



7 



More precisely, let 









/ 7o 


7i 




7fc-i^ 






— 1, /ijfc = det 


7i 


72 




7fc 


, k> 1, 








\7fc-i 


7fe 




l2k-2 J 








/ 7o 7i ■ 




lk-2 


7fc \ 




h k 


= det 


7i 72 ■ 




7fc-i 


7fc+i 


jfe > 2. 






\7fc-i 7fc • 




72fc-3 







Then 



hk-i hk+i 



1/2 



7 
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Since 



({H-zi)- i 5 X ,8 l )=-- l ttm 



where H is given by (|7.3j) . the function tan f~) /ios the following Taylor 
expansion at infinity 



tan 



— 2 2n (2 2 " - 1) 5 ; 



2 1 1 



71=1 



(2n)! z 2 ™" 1 ' 



where {B2 n } are Bernoulli numbers [21], B2 = 1/6, B^ = —1/30, B§ = 
1/42, . . we get that 

Z 2( fc -l) 22fc(2 2 fc _ 1)gafc 

72fc = 7777^ , 72*-i = 0, fc = 1, 2, . . . . 

(2*)! 

The formula for {a^} allows to get some relations between Bernoulli 
numbers. 

Theorem 7.5. Let J t be a Jacobi matrix of the form 
1 



( 1 ~hi 00 



l n I 



-A? 



1-3 V3-5 

* n I 







'3-5 V5-7 



V 



i ' 

y'(2n-3)(2n-l) v /2n-l)(2n+l) 

7 



where t is a complex number and I is a positive number. Then the 
corresponding Jacobi operator for all t but two t = ±il has a complete 
system of eigenvectors. 

Proof. The matrix J t has the representation 

J t = H + t(-,5 1 )5 1 , 

where H is given by (|7.3|) . Suppose that Imt 7^ 0. Then the non-real 
eigenvalues of J t are solutions of the equation 

l + tdH-ziy^udx) =0. 

Because 

((H-zI)- 1 6 u 6 1 ) = ~taa 
we get the equation 

. l\ I 

tan 
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This equation has the following solutions: 



21 ( arg 



t + il 
t-il 



+ 2irk + i In 



t + il 



t - il 



arg 



t + il 
t-il 



+ 2vk\ + In 2 



t + il 



t-il 



, keZ, t^ ±il. 
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Hence 



ImZk = Imt. 
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Since in this case the eigenvalues of Jt have the algebraic multiplicities 
one, according to the Livsic Theorem 12 .31 the operator J t has complete 
system of eigenvectors. 
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If Imt = 0, then the operator J t is self-adjoint and compact, its 
eigenvalues are real numbers 

Zk = 77— -,\ , k e Z, 



and the corresponding eigenvectors form a complete system. □ 
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